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ABSTRACT 

Aims. We constrain the dust distribution and its properties (temperature, emissivity) in inner proto-planetary disks 

Methods. We performed sub-arcsecond high-sensitivity interferometric observations of the thermal dust emission at 1 .4 mm and 2.8 mm in the 
disks surrounding LkCa 15 and MWC480, with the new 750 m baselines of the IRAM PdBI array. This provides a linear resolution of ~ 60 
AU at the distance of Taurus. 

Results. We report the existence of a cavity of ~ 50 AU radius in the inner disk of LkCa 15. Whereas LkCa 15 emission is optically thin, the 
optically thick core of MWC480 is resolved at 1.4 mm with a radius of ~ 35 AU, constraining the dust temperature. In MWC480, the dust 
emission is coming from a colder layer than the CO emission, most likely the disk mid-plane. 

Conclusions. These observations provide direct evidence of an inner cavity around LkCa 15. Such a cavity most probably results from the tidal 
disturbance created by a low-mass companion or large planet at ~ 30 AU from the star. These results suggest that planetary system formation 
is already at work in LkCa 15. They also indicate that the classical steady-state viscous disk model is too simplistic a description of the inner 
50 AU of "proto-planetary" disks and that the disk evolution is coupled to the planet formation process. The MWC 480 results indicate that a 
proper estimate of the dust temperature and size of the optically thick core are essential for determining the dust emissivity index /?. 
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1. Introduction 

With current instrumentation, proto-planets located at the dis- 
tance of Taurus cannot be observed directly. The best evidence 
of planet formation in proto-planetary disks remains the detec- 
tion of the tidal gaps. Because of the dependence of the dust 
temperature on the distance from the star, these gaps should 
leave signatures in the infrared SED of the objects, in the form 
of a deficit of em i ssion at wavelengths between ~ 1 and 30 
/mi. Koerner et al. ( 19931) suggest that the deficit of emission at 
10 /im in GM Aur was indeed the result of such a gap. Many 
recent studies of the spectral energy distribution (SED) in the 
near-infrared (NIR) and mid-infrared (MIR) coming from data 
obtained by various telescopes, in particular from the Spitzer 
satellite, have revealed NIR-MIR emission dips that can be in- 
terpreted by tru ncated disks with i nner radii of ~ 3-5 AU 
up to ~ 25 AU (ICalvet et al.ll2005l) . Some direct evidence of 
inner cavities has been also obtained in the scattered light im- 
ages of the disks of moderat e opacity associated w ith Herbig 
Ae stars such as HD141569 dAugereau et al.lll999l) . More re- 
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cently, a large inner radius (70 - 100 AU) was measured in 
the mm domain for the d ense spiral-like structure found in AB 
Auriga jPietu et al.ll2005l) . 

However, models of planet formation predict gaps with 
density contrasts of only ~ 10-100 for a Jupiter-like planet (e.g 
Crid a et al. I l2006l) . Due to the very high dust opacity observed 
in the IR (at ~ 1/rni, assuming standard dust g rain properties, 
the o pacity can be as high as 10 4 - 10 6 , see iD'Alessio et al 



19981) . the IR SED does not allow astronomers to get detailed 
knowledge on the moderate de nsity contrasts. Th i s can b e qual- 
itatively understood in Fig. 5 in lDullemond et alJ (1200 II) . where 
a change in the slope p of 1 on the dust surface density does 
not significantly affect the SED (by less than ~ 20 %), which is 
therefore difficult to disentangle from the other physical effects 
on real data. 

As a consequence, only gaps corresponding to a significant 
decrease in the surface density in the inner disk, i.e. a transi- 
tion in opacity between thick and thin regimes, can be detected 
in the modelling of the IR SEDs. Only resolved images ob- 
tained with a tracer of moderate opacity can provide reliable in- 
sights into the geometrical structure and density contrasts when 
the inner cavities are not significantly empty. The dust opacity 
of proto-planetary disks in the millimeter/sub-millimeter range 
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Fig. 1. Top row: Results for MWC 480. From left to right: 2.8 mm continuum image. The angular resolution is 1.09 x 0.74 at 
PA 28°, and the contour spacing is 1 mJy/beam (0.13 K, about 3.5 <x). The 1.4 mm continuum image: the angular resolution is 
0.73 x 0.53 at PA 33°, and the contour spacing is 5 mJy/beam (0.33 K, 2.5cr). The 1 .4 mm residual from the best fit with the same 
contours. Right: 1.4 mm real part of the visibility vs baseline lengths. The points and error bars represent the measured values, 
and the curve indicates the visibilities of the best-fit disk model. The visibilities have been corrected from the disk orientation 
and inclination by compressing the V value by cos(/) prior to circular averaging. Bottom row: Results for LkCa 15. The 2.8 mm 
image is obtained by considering baselines longer than 200 m (but no flux was lost in this process). The resolution is 1.06 x 0.61 
at PA 31°, and the contour spacing is 0.75 mJy/beam (0.12 K, 2 &). The 1.4 mm resolution is 0.66 x 0.39 at PA 55°, and the 
contour spacing is 1.5 mJy/beam (or 0.15 K, 1.7 cr). The crosses indicate the position, orientation, and aspect ratio of the disks. 



is well adapted to imaging such "young" cavities, provided 
enough angular resolution and sensitivity is available. 

In this paper, we present new PdBI observations at an angu- 
lar resolution of 0.35" - 0.5" for the Keplerian disks surround- 
ing two bona-fide single PMS stars LkCa 15 and MWC 480. 
LkCa 15 is a classical TTauri star of spectral type K5, mass 
~ 1 M , and age ~ 3 - 5 Myr and MWC 480 is a Herbig Ae 
star o f spectral type A4 and age 
20001) . 



5-7 Myr dSimon et al 



2. Observations and results 

The observations were performed from winter 2001 to winter 
2005. We observed simultaneously at 110 and 220 GHz (2.8 
and 1.4 mm). Typical system temperatures range from 200 to 
400 K at 220 GHz. We mainly used 4 configurations: the stan- 
dard BCD configuration of the IRAM array, with (projected) 
baselines ranging from 15 m to 320 m, and the new A + con- 
figuration, which provides baselines extending up to ~750 m. 
Some additional data were obtained on MW C 480 in A config - 
uration, using time sharing with AB Aur (see Piet u et al . 2005). 
The data were obtained under good weather conditions: the rms 
phase did not exceed ~ 70° at 1 .4 mm even on the longest base- 
lines. The two sources were observed in track-sharing mode, 
alternating every 10 minutes from one source to the other. For 
the BCD configurations, we calibrated the data using the stan- 



dard method, with MWC 349 as a reference for flux calibration. 
On the longest baselines, MWC 349 is partially resolved at 1 .4 
mm, and a different method was required for the flux calibra- 
tion. We first based our flux scale on an extrapolation in time 
of the flux densities of the two phase calibrators, 0415+379 
(i.e. 3C111) and 0528+134. However, since quasars have un- 
predictable flares, we also compared the measured flux density 
of MWC 480 obtained from the BCD data to that obtained in 
the A + data on the overlapping baseline range (from 100 to 
400 m). MWC 480 is bright and compact enough to allow such 
a direct comparison, to a precision of ~ 10 %. The two methods 
give the same results within 5 %. The simultaneous observation 
and common calibration ensure that any morphological or in- 
tensity difference observed between the two sources is real and 
not due to instrumental artifacts. 

At 220 GHz, these observations provide an angular reso- 
lution about 0.7 x 0.4" using natural weighting. With uniform 
weighting, the longest baselines (500 kA) allow us to reach an 
angular resolution of about 0.5x0.3" at PA ~ 45°. At 110 GHz, 
the resolution is a factor lower. Based on the integration time, 
system noise, and measured efficiencies, the expected (thermal) 
noise level was 0.7 mJy/beam at 220 GHz. However, the dy- 
namic range is limited by phase noise. This results in an effec- 
tive noise of 0.9 mJy/beam for LkCa 15 and 2.0 mJy/beam for 
MWC 480. At 110 GHz, the noise is 0.3 mJy/beam, so essen- 
tially thermal. 
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Table 1. Comparison of the inclinations and orientations of the 
LkC a 15 and M WC 480 disks derived from CO isotopologues 
(see Pietu et al. 2006) and mm continuum emission. 





from CO from 2.8 and 1.4 mm 




isotopologues 


emission 




LkCa 15 




PAD 


150 ± 1 


151 ± 3 


iC) 


52 ± 1 


49 ±3 




MWC 480 




PAD 


58 + 1 


61 ±8 


«o 


36+1 


31 ±5 



Table 2. Parameters of the MWC 480 disk derived from the 
model fits. 



Fixed input parameters 




tf,n, (AU) [3f 


[3] 


q [0] 


[0.5] 


Derived parameters from 2.8 and 1.4 mm dust emission 


r 100 (K) 24.5 ± 3 


9.8 ± 0.4 


R out (AU) 190 ± 15 


185 ± 15 


£ioo(g.cnr 2 ) 3.1+0.2 


17.1 ± 1.5 


p 2.5+0.1 


1.6 + 0.1 


P 1.26 + 0.05 


1.51+0.06 


M disk (M Q ) 0.24 


0.30 


M(fl>35AU)(M Q ) 0.04 


0.19 



The error bars are the lcr formal errors (effective noise, see Sect.2) 
from the fit. R int has been fixed to 3 AU. (a) Square brackets indicate 
fixed parameters. 

Figure Q] presents the images obtained at 2.8 mm and 
1.4mm on MWC 480 and LkCa 15. While the disk of 
MWC 480 appears centrally peaked, the disk around LkCa 15 
does not. The 1.4 mm imag e of LkCa 15 is rather similar to 
the first image of GG Tau dDutrev et alJll994 , see their Fig. 1) 
where the inner radius of the circumbinary dust ring was just 
resolved by the interferometer. The total flux densities at 2.8 
and 1.4 mm are 17.0 + 0.8 mJy and 140 ± 3 mJy, respectively, 
for LkCa 15, corresponding to an apparent spectral index of 
a = 3.04 ± 0.07. For MWC 480, S(2.8) = 35.2 + 0.8 mJy 
and S(1.4) = 235 + 4 mJy, leading to a = 2.14 + 0.04. The 
calibration uncertainties result in an additional error of ±0.15 
on the absolute values of a, but the relative values are not af- 
fected, since the data were obtained simultaneously. Figure Q] 
also shows the residual images after subtraction of the best- 
fit model (see Tables |2][3]). The rightmost panels of Fig[TJ dis- 
play the real part of circular average of the calibrated visibilities 
obtained after deprojection, superimposed with the visibilities 
from the best models (curves). LkCa 15 is well resolved with a 
first null at 130 kA (and probably a second null at ~ 400 kA). 

3. Data analysis 

To better quantify the obvious morphologic differences be- 
tween both sources, the 2.8 mm and 1.4 mm data were ana- 
lyzed in the Fouri er Plane using our disk fitting method (see 



Table 3. Parameters of the LkCa 15 disk derived from the 
model fits. 

Temperature from CO isotopologues 

r 100 (K) 22+1 
_q 0.37 ± 0.02 

Derived parameters from 2.8 and 1.4 mm dust emission 

R mt (AU) 46 ± 3 [5]" 

R oui (AU) 177 ± 12 135 ± 3 

S 100 (g.cm- 2 ) 3.1+0.4 4.1+0.6 

p 1.7 + 0.3 -0.5 + 0.2 

P 1.2 + 0.1 1.1+0.1 

M disk (M G ) 0.029 0.025 

The error bars are the lcr formal errors (effective noise, see Sect.2) 
from the fit. Two solutions are presented: one with a free inner radius, 
and one with a fixed inner radius of 5 AU. The dust emission is too 
optically thin to allow us an independent fit of the tempera t ure, w hich 
is adopted from the CO analysis presented in IPietu et al] ( 120061) . (a) 
Square brackets indicate fixed parameters. 



law distribution of the temperature T(r) = Tioo(r/lOO AU) 9 , 
surface densit y £(r) = £ioo(r /TOO AU)~ P . Follo wing 



Dutr ev et alj |2006, and references therein). We modelled an 



inclined disk with truncated inner and outer radius and power 



Beckwith et"all dl990l) . we used the iHildebrandl d 1983b pre- 
scription for the dust opacity, k(v) = O.^v/lOOOGHzy 3 cm 2 /g 
(per gram of dust+gas). We simultaneously fit the 2.8 and 
1.4 mm UV data, thereby constraining (3 directly. As shown 
in Table Q] for both sources, the inclinations and orienta- 
tion derived from the continuum data agree with those deter- 

19 II — 

mined from an analysis of the CO (Simon et al. 2000) and 

I3 CO observations (for details, see lPietu et al.ll20 06). We used 

the more accurate CO-based inclinations in the analysis. For 

MWC 480 the temperature can be constrained by the dust ob- 

serv ations. For LkCa 15, we used the temperature determined 

by IPietu et all d2006l) from CO lines. We checked that the 

choice of the scale height does not influence the results at all: 

we used h(r) - 16.5(r/100 AU) L25 AU for both sources. 

For continuum data, the brightness distribution scales in the 
optically thick case as the temperature T, and its exponent q, 
while in the optically thin case, it scales as the product "LT and 
the sum p + q. Since in general the dust opacity decreases with 
radius, one can measure both the surface density and the tem- 
perature with angular resolution that is high enough. 

Table E] shows the best fits for MWC 480. The large value 
of the peak surface brightness (7 K in FigJTJ indicates a signifi- 
cant opacity at both wavelengths in the inner ~ 35 AU. Indeed, 
optical depth reaches unity between r — 31 AU (model 1) and 
38 AU (model 2) at 220 GHz, and between 20 and 23 AU at 1 10 
GHz. We thus attempted to treat the dust temperature as a free 
parameter. The measurements directly constrain the tempera- 
ture around 20 - 30 AU to be ~ 20 K. However, the optically 
thick core is not sufficiently resolved to determine the temper- 
ature exponent q. Table [2] shows two solutions. The data are 
compatible with a Tioo = 24 K, q = or with steeper laws, 
^ioo - 10 K, q — 0.5, but not with higher temperatures. The 
temperature derived from 12 CO (T^o = 45 K and q - 0.65) is 
much too high to apply to the dust. All solutions intersect near 
r = 17 AU with T = 24 K. The lack of knowledge on the tem- 
perature law results in a significant uncertainty in the surface 



4 



Pietu et al.,: Resolving the LkCa 15 and MWC480 inner disks 



density law, and thus on the total disk mass. With these opaci- 
ties, the optically thick inner core contributes quite significantly 
to the total flux, and the dust emissivity index (3 — 1 .4 + 0.15 is 
thus substantially higher than a -2. 

For LkCa 15, the surface brightness of the mm emission 
is low (< 3 K), indicating that the dust is essentially optically 
thin. It is thus impossible to constrain the dust temperature with 
this data. On t he o ther hand, an analysis of I2 CO and I3 CO 
by iPietu et al.l (120061) provides a unique value for the temper- 
ature, Tioo =22+1 and q = 0.37 + 0.02. We adopted these 
values. With a small inner radius (5 AU, but the result do not 
change for lower values), the best-fit solution implies p -0.5, 
i.e. a surface density increasing with radius. This is physically 
implausible. Treating the inner radius as a free parameter re- 
sults in a better fit at the 4 <j level and a more physical value 
for p (see Table |3J. Other combinations of inner, outer radius 
and density exponent are possible: for example, a solution with 
p - 3, Rint - 65 AU and R out > 600 AU (in agreement with the 
disk extent in CO) is only 2-3 cr above the best-fit solution for 
LkCa 15. This may indicate that LkCa 15 is not devoid of dust 
beyond ~ 170 AU, but surrounded by a more tenuous disk. 
Despite the correlation between p and R m , the LkCa 15 disk 
exhibits a clear decrease in surface density towards its center, 
with the best solution indicating a 50 AU radius central hole. 

In conclusion, while MWC 480 presents a "classical" dust 
distribution, LkCa 15 indicates that dust disks can present 
brightness distributions versus radius that are more sophisti- 
cated than simple power laws, and suggests the existence of a 
large inner cavity in an otherwise "standard" dust disk. 



The value of ft = 1.4 + 0.15 indicates moderate grain 
grow th. This i s sim ilar to the value obtained for AB Aur 
by IPietu et aT (2005). Using the coronagraphic mode o f the 
HST/NICMOS camera at 1.6//m, lAugereau et"al] bOOlb ob- 
served the disk of MWC 480, but did not detect its scattered 
light. Their upper limit on the column density of scattering ma- 
terial is ~ 10~ 4 g.cirT 2 at 100 AU (of dust only), several orders 
of magnitude below our measured value. Our results suggest 
that a significant amount of dust is still "hidden" to the scattered 
light regime, most likely because the dust disk is geometrically 
thin and the starlight partially masked by the dust orbiting close 
to the star. This can occur if the dust has significantly settled 
on the disk-mid plane. This explanation a grees with the clas- 



sification of MWC 480 as a g roup II star dMeeus et all 12001 



Acke & van den Anckerl2004 ). The constraint on the dust tem- 



perature indicates that the dust is colder than the zone traced 
by the CO isotopologues, in further agreement with this hy- 
pothesis. Although indirect, these arguments suggest that the 
dust has settled towards the disk plane. Table [2] indicates the 
MWC 480 disk is quite massive. However, this is largely the 
result of extrapolation of the density and/or temperature laws, 
as indicated by the total mass beyond 35 AU in Table[2] 

MWC 480 illustrates the need to resolve the inner opti- 
cally thick core and to measure the dust temperature to prop- 
erly characterize t he dust emissivity. Lower-resolution data 
(IDutrey et all 19961) typically tend to find values of (3 1, and 
our new results indicate such values could be significantly bi- 
ased by contamination from optically thick dust at low temper- 
ature. 



4. Discussion 



4.2. The inner disk of LkCa 15 



4.1. MWC 480 

MWC 480 was recen tly observed at 1.4 mm by 
d2006l hereafter HLM) using the BIMA 



Hamidouche et al 



array, with slightly higher resolution, but lower sensitivity. 
HLM mention a similar source orientation. They use a temper- 
ature derived from the IR SED, T\oq = 17 K, and q = 0.62. 
However, because of the IR opacity, the IR-constrained 
temperature deduced from an isothermal vertical disk model 
naturally overestimates the mid-plane temperature to which 



the mm emission is sensitive (see also iD'Alessio et al.M 1998 



for details). Indeed, the temperature used by HLM is much 
too high to explain our observations, and using it results in 
a best-fit solution that is llcr worse than those presented in 
Table [2] HLM also derive a relatively shallow surface density 
distribution, p = 0.5 - 1.0, as a result of the high assumed 
temperature. The HLM temperature law requires a low opacity 
to reproduce the surface brightness (9 K in their Fig.l). As a 
consequence, the shallow inner-surface-brightness distribution 
is represented by a flat surface density, while it is in reality due 
to contribution by the optically thick core. Also, p is biased 
towards low values by the extended emission present in the 
HLM image, which is not confirmed in ours. This structure 
may be due to phase errors. Note that in general, phase noise 
will result in a flattening of the brightness distribution, since it 
scatters signal (much like an effective seeing). 



Until now, observations from the NIR to the mm range have 
shown that LkCa 15 is a bona-fide Classical TT auri star sur- 
rounded by a large C O disk of a few .01 M ( Simon et al 



l2000b lOi et alJl^Oolh . iLeinert et alJ (1 19931) observed LkCa 15 
and classified it as a single star, since they found no companion 
down to a separation of ~ 0.13" (20 AU). 



Bergin et al. (2004) have observed the H2 UV emission and 



the NIR-MIR SED of LkCa 15. They conclude that the disk has 
an inner radius (puffed-up inner rim) located at ~ 3 AU. From 
their disk modelling, they estimate a column density of dust 
^dust - 0.2g.cirT 2 at 1 AU, assuming a standard accretion disk 
model with a = 0.01 and T = 100 K. Extrapolated to 45 AU, 
this translates into 0.07 g.cirT 2 of dust+gas (with a gas-to-dust 
ratio of 100). This is 200 times lower than the surface density 
from Table [2] confirming the interpretation of the mm images 
of LkCa 15 as the result of a cavity of radius ~ 45 AU nearly 
devoid of dust. If we fit a p = 1 surface density law into the 
LkCa 15 cavity, we find an upper limi t of Zdust - 0-12 g. cm' 2 



at 1 AU, consistent with the results of Berg in et al.l ((2004). 

In the current knowledge of inner-disk destruction pro- 
cesses, among the mechanisms (such as planetary formation, 
photo-evaporation, radiation pressure, etc. . . ) that can explain 
this distribution, planetary formation or a low mass companion 
remain plausible since the LkCa 15 disk is fairly massive. The 
large disk density makes photo-evaporation and radiation pres- 
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sure unlikely to be effective. lAlexander et al. I d2006l) showed 
that when photo-evaporation starts propagating beyond 20 - 30 
AU, the surface density near 50 AU is 0.03 g.cirT 2 , within 
a factor of a few to account for different UV fluxes (see their 
Fig.l). This is 500 times lower than derived from our best fit 
model in Table [3] 

In a viscous disk, a gap can be opened by a planet as soon 
as the Hill's radius is larger than the disk height. Following 



Take uchi et al. (1996), the half-width of a gap created by a 



proto-planet or a low-mass companion can be approximated 
by: w = 1.3aA'/ 3 where a is the semi-major axis of the orbit. 
The strength ratio of tidal to viscous effects A is : 



A = (M p /M t y 



1 



3a(h(r)/r) 2 



For a 30 AU and assuming a ~ 0.01, h(r)/r ^ 0.1, M t 
1 M and the mass of the planet M p m 0.005 M , we find A 
0.083 and w - 17 AU. Hence a ~ 5 - 10 Jupiter mass planet 
orbiting at 30 AU would be sufficient to evacuate the inner 50 
AUoftheLkCal5disk. 

Getting an upper limit on the or biting object mass is more 
difficult. The separation limit from iLeinert et al.l ( 1993 ) is 20 
AU, but given the inclination of the LkCa 15 disk, an ob- 
ject orbiting at 30 AU would spend most of its time at pro- 
jected distances about 20 AU. The best limit on th e companion 



mass comes from the kinematic determination of Simon et al 



(2000), who indicate a total mass for the system of 1 .0+0. 1 Mq. 
The spectral type of LkCa 15 makes it unlikely to be a star of 
less than 0.8 M , leaving at most 0.2 M for a companion. A 
coeval star of this mass would have a K magnitude of about 
10.5. 

5. Summary 

We used the new 750 m baselines of the IRAM array to image 
the dust disks orbiting MWC480 and LkCa 15 at an angular 
resolution of ~ 0.4". 

- The dust disk of MWC 480 is centrally peaked with an opti- 
cally thick core of radius ~ 35 AU at 1 .4 mm. Since the core 
is resolved, we can estimate its temperature, which is sig- 
nificantly lower than that of the CO layers, and measure the 
dust emissivity index ft. The low dust temperature suggests 
that grains emitting at mm wavelengths tend to be located 
around the colder disk mid-plane, in agreement with sedi- 
mentation. The index f3 is ~ 1 .4, suggesting moderate grain 
growth, like in AB Aur. 

- To be compatible with existing NIR images of MWC 480, 
the dust must be confined to a geometrically thin disk and 
mostly hidden from the starlight. 

- The LkCa 15 disk reveals a large (~50 AU radius) cavity, 
which is not completely devoid of dust. The most proba- 
ble explanation for such a wide cavity is the existence of a 
massive planet (> 0.005 Mo) or of a low mass companion 
(< 0.2 Mo) orbiting around ~ 30 AU from the star. 

Long baselines on large mm/submm arrays are now opening 
the planetary formation regions of nearby young disks to in- 
vestigations. However, as for AB Aur, figuring out whether the 



companion to LkCa 15 is a proto-planet or a low-mass star re- 
quires deeper IR images, or the advent of ALMA, which may 
resolve the structure of the cavity. 
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